Bounds for the second and the third derivatives of the electron density at the nucleus 
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Lower bound for p"(0), the second derivative of the spherically averaged atomic electronic density 
at the nucleus, and upper bound for p"'(0), the third derivative, are obtained respectively. It is 
shown that, for the ground state, p"(0) > ?§-Z 2 p(0) and p"'(0) < -^Z 3 p(0) where Z is the charge 
of the nucleus, and p(0) is the electron density at the nucleus. Tighter bounds for p"(0) and p"'(0) 
are also reported which are valid for both the ground state and excited states. Explicit illustration 
with the example of one-electron atomic ions is given. 

PACS numbers: 



I. INTRODUCTION 

Rigorous knowledge of the electron density is rather 
valuable in understanding the electronic structure of 
atoms, molecules, and solids. It is also helpful in guiding 
the construction of the approximations for the exchange- 
correlation energy functional and the correspondingpo- 
tential in the approach of density functional theory [l[ to 
the problems of inhomogeneous electron systems. One of 
the well known exact results is the so-called Kato theo- 
rem Q which states that, near any nucleus, 



p'{r)\ r = = -2Zp{0). 



(1) 



Here Z is the charge of the nucleus which is taken as 
the origin; f(r) means the spherical average of function 
/(r); primes denote the derivatives with respect to r in 
this paper. 

In recent work Q , by investigating the behavior of the 
wavefunctions of the interacting Schrodinger system and 
the corresponding noninteracting Kohn-Sham system in 
the vicinity of the nucleus, we established relations be- 
tween the second and the third derivatives of the spheri- 
cally averaged density at the nucleus. They could be un- 
derstood as extentions of the cusp condition of Eq. (Q]) 
to higher orders of derivatives. In this paper, we derive 
rigorously a lower bound for the second derivative and 
an upper bound for the third derivative, respectively, of 
the spherically averaged density at the nucleus for the 
ground state. The bounds are given as follows: 



and 



p"(0) > yZ 2 p(0), 



p"'(0) < -yZ 3 p(0). 



(2) 



(3) 



The results hold in general whether the system is an 
atom, molecule or a solid. Tighter bounds for p"(0) and 



p"'(0) are also reported (see Eqs. ([23]) and ([24]) in Sec. 
Ill) which are valid for the excited states as well. 

In Sec. II, we discuss the near nucleus behavior of 
the wavefunction of the Schrodinger equation and the 
density. The derivation for the bounds is presented in 
Sec. III. In Sec. IV explicit illustration for the case of 
one-electron atomic ions is given. Summarizing remarks 
are made in Sec. V. 



II. NEAR NUCLEUS BEHAVIOR OF THE 
WAVEFUNCTION AND THE DENSITY 

The Schrodinger equation for TV-electrons in an exter- 
nal potential v(r) arising from their interaction with nu- 
clei is (in a.u.) 



= (f + V + U)^ = E^, 



(4) 



where t = -|EiV?, V = E< & = 
lEi^j \n-r- \ ' ^ * ue wavefunction, and E the energy. 
Following Ref. Q, we write, for limiting small r, the 
many-body wavefunction as 

*(r, X) = *(0, X) + a(X)r + b(X)r 2 + c(X)r 3 + . . . 



l=— 1 
2 

[b 2m (X)r 2 + . . .]Y 2m (r) 



m=-2 



(5) 



where f = r/r, X denotes s,r2S2, ■ ■ ■ ,?nsn, and Yi m (f) 
are the spherical harmonics. For r — > 0, the Schrodinger 
equation |gj) can be rewritten as [3j 



rn — — 1 



H^(X)*(r,X) = E*(r,X), 



where 
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N 

9m ( X ) = ^^-^y 1 * ro (f i ), 



3 ^ r- 

i=2 * 



(6) 



(7) 



2 



and 



fff_- 1 1 (X)= ^[-i v 2 +u(r ,) + I] 

i=2 1 

Substituting the expression of Eq. ([5|) into Eq. (J6j> and 
equating the coefficients of the terms of r" 1 , r°Yi m (f), 
r , and r , respectively, one has 

a(X) + Z*(0,X) = 0, (9) 
26 lm (X) + Za lm (X) =0, (10) 
36(X) - Z 2 *(0,X) = [iJ^J^X) - B]*(0,X),(11) 
6c(X) + Z6(X) = Z[E - i^^X^^X). (12) 

We note that the behavior of the wavefunction in the 
vicinity of the nucleus has been extensively investigated 
In fact, Eqs. (9) and (10) have been obtained 
previously. On the other side, Eqs. (11) and (12), which 
shall play a key role in the derivation in Sec. Ill, have 
not been previously reported in this separate form. 
The density is defined as 



p(r)=N J dX|*(r,X)| 5 



(13) 



where J dX denotes Yl s I °^ X2 ■ ■ - dX-N- With the rela- 
tions in Eqs. (9) and (10), it can be shown that 

p"(Q) = 2Z 2 p(0) + AN J dXRe[»f*{0,X)b(X)} 

r 1 i 

+2N JdXj2 ^|ai m (X)| 2 (14) 



and 



p'"(0) = 12N J dXRe[V*(0,X)(c(X) - Zb{X))] 

r 1 1 
-6ZN J dXj2 ^l«i™( X )! 2 , (15) 



m= — 1 



(see also Ref. @.) This completes the discussion of the 
behavior of the wavefunction and the density near the 
nucleus. 



III. DERIVATION OF EQS. (HJ AND ^ 

We observe that Eqs. (11), (12) and Eq. (TJ3]) lead to 
N J dX Re[$*{0,X)b{X)} 

> \{Z 2 -E + E N z zl )p{Q), (16) 



N J dX fte[#*(0,X)c(X) 



< -Z{-Z* + AE-AE»lt fi )p{Q), (17) 



where E z _l Q is the ground state energy of the Hamil- 

tonian H Z Z{{X). Inequalities JT6J) and (HH) are critical 
equations in the following derivation. Comparing them 
with Eq. (TT4"|) and Eq. (TT5)) . respectively, one has 

p"(o) > ^V(o) + \{E N z :i Q - e) p (o) 
r 1 i 

+2N dX ^l«i™( X )| 2 , (18) 

ra— — 1 

and 

p"'(Q) < -y Z*p(0) - - E)p(0) 

r 1 i 

-6ZN dX ^|ai m (X)| 2 . (19) 
J m=-l 

In passing, we mention that, at the nucleus, the kinetic 
energy density, which is defined as 

t(r) = ±N J rfX V **(r,X). V *(r,X), (20) 

has been shown in Ref. Q as 

t(0) - \z 2 p{Q) = NjdXj2 ^\*UX)\ 2 . (21) 



Equation (f2~Tj) indicates that 

1 



t(0) > z Z 



(22) 



The existence of the term on the right hand side of Eq. 
0} had not been recognized before Q. 
From Eqs. (JT5J) and (fTO)) , one obtains 

p"(0) > fz 2 p(0) + \{E N z zl, - E)p(0), (23) 



and 



p"'( ) < ~Z 3 p(0) y Z^r^o - ^)P(O). (24) 

Up to this point, all the calculations are in fact not re- 
stricted to the ground state. The bounds shown in Eqs. 
(|23|) and (|24|) are valid for both the ground state and ex- 
cited states. For the ground state we further obtain Eqs. 
© and © from Eqs. fl23J| and 
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IV. ILLUSTRATION FOR ONE-ELECTRON 
ATOMIC IONS 



The inequalities of Eqs. ([25)) and ([26)) are satisfied. No- 
tice that in this case a\ m ^ in Eqs. (fl8|) and (|19p . 



For one-electron atomic ions (including hydrogen 
atom), Eqs. and ((Ml) become 



and 



p"(0)>-(5Z 2 -2E)p(0), 



p'"(0) < --^(7Z 2 - 10E)p(0). 



(25) 



(26) 



The wavefunctions are exactly known for both the ground 
state and excited states: *(r) = R„i(r)Yi m (r), with @ 



Rni(r) = - 



,3+21 



(n-l-l)\ 
2n[(n + l)\} 3 



1/2 



r/2 



r l L%+t(a n r), 



(27) 



where a n = 2Z/n, and L^^(a n r) are the Laguerre poly- 
nomials. It is easy to see that for 1^1, a\ m = in Eqs. 
(fl"8f and (19]), and Eqs. (23]) and d24j) in fact become 
equalities: 



and 



p"(0) = -(5Z 2 -2£)p(0), 



p"'(0) = -^Z(7Z 2 10E) P (0). 



(28) 



(29) 



We next give an explicit illustration of the results with 
the exactly known wavefunctions. Obviously for I > 2 
Eqs. (HHJ) and (JUJ) are trivially true since all p(0), p"(0), 
and p"'(0) are zero. For I = 0, one can obtain from Eq. 
(123 



Z 3 

p(r) = — o 



1 



1 



2Zr+ -Z 2 (5 + -^)r 



1 73/7 , ^ \„3 



-z%7 + ^y 



(30) 



Equations ([28)1 and ([29)) are hence confirmed with the 
help of the fact that E = —Z 2 /2n 2 . For / = 1, one has 



p(r) 



1 

9tt \ n 



(n 2 -\)r 2 {\- Zr + ...). 



(31) 



V. SUMMARY 

In summary, in this paper, we have established the 
lower bound for p"(0), the second derivative of the spher- 
ically averaged electron density at the nucleus, in Eq. 
(2), and the upper bound for p"'(0), the third deriva- 
tive, in Eq. (3). Tighter bounds for p"(0) and p"'(0), 
valid for both the ground state and excited states, are 
also reported in Eqs. ([23)1 and ([24| . These results add 
some to our rigorous information of the electron density, 
which might be valuable in the calculation of the elec- 
tronic structure of atomic systems. 

Whether the opposite bounds exist for p"(0) and p"'(0) 
remains an interesting question. Not unrelated to this 
question is the fact that both upper and lower bounds 
have been extensively explored for p(0) (and, according 
to Eq. ©, equally for p'(0)) 0- 

Note added. During the preparation of this 
manuscript, the author became aware of related work 
Q. The bounds shown in Eqs. ([2]) and ([3]) in this paper 
are indeed the same as those in Eqs. (1.25) and (1.32) in 
Ref. (up to a different form of Hamiltonians). There 
are, however, minor differences between inequalities ()23p 
and ([24)) of the present paper and the first inequalities of 
Eqs. (1.25) and (1.32) in Ref. 0. In fact, E%zl Q ~E can 



be rewritten as E^_l Q - 



E 



N-l 



N-l 



■ p, where p, = E — E z 
is the minus the ionization energy and denoted as — e in 



Eq. (1.20) in Ref. 



Since E%zt fi 



E 



N-l 
Zfi 



> 0, it is 



easy to see that the bounds in Eqs. ([2"3")l and ([M)l are 
tighter than those given by the first inequalities of Eqs. 
(1.25) and (1.32). The differences were also recognized 
and analyzed in Remark 1.4 in Ref. 0, though equations 
like inequalities ([2"3"| and ([2"4")l were not explicitly given 
there. 

The author is indebted to Profs. S. Fournais, 
M. Hoffmann-Ostenhof, T. Hoffmann-Ostcnhof, and T. 
0stergaard S0rensen for confirming that the bounds 
shown in Eqs. ^ and ^ are the same as those in Eqs. 
(1.25) and (1-32) in Ref. [§]. Enlightening comments on 
this paper from these colleagues are also gratefully ac- 
knowledged. This work was supported by the National 
Science Foundation of China under Grant No. 10474001. 
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